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ABSTRACT 



We calculate the statistical distribution of observed afterglow flux, in some 
fixed observed frequency band, and at some fixed observer time after the 
explosion (t bs) in t wo models - one where the explosion takes place in a uniform 
density medium and the other where the surrounding medium has a power-law 
stratification such as is expected for a stellar wind. For photon energies greater 
than about 500 electron-volt and t obs > 10 3 sec the afterglow flux distribution 
functions for the uniform ISM and the wind models are nearly identical. We 
compare the width of the theoretical distribution with the observed x-ray 
afterglow flux and find that the FWHM of the distribution for energy in 
explosion and the fractional energy in electrons (e e ) are each less than about one 
order of magnitude and the FWHM for the electron energy index is 0.6 or less. 



Subject headings: gamma rays: bursts - gamma-rays: theory 
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1. Introduction 

The improvement in the determination of the angular position of gamma-ray bursts 
(GRBs) by the Dutch- Italian satellite Beppo/SAX has led to the discovery of extended 
emission in lower energy photons lasting for days to months, which has revolutionized our 
understanding of GRBs (cf. Costa et al. 1997, van Paradijs et al. 1997, Bond 1997, Frail et 
al. 1997). The afterglow emission was predicted prior to their actual discovery by a number 
of authors (Paczyhski & Rhoads, 1993; Meszaros & Rees, 1993; Katz, 1994; Meszaros & 
Rees 1997) based on the calculation of synchrotron emission in a relativistic external shock. 
The afterglow observations have been found to be in good agreement with these theoretical 
predictions (cf. Sari, 1997; Vietri, 1997a; Waxman, 1997; Wijers et al. 1997). 

The medium surrounding the exploding object offers some clue as to the nature of 
the explosion. Vietri (1997b), and Chevalier and Li (1999a-b) in two very nice recent 
papers, have offered evidence that some GRB afterglow light curves are best explained 
by a stratified circumstellar medium which suggests the death of a massive object as the 
underlying mechanism for gamma-ray burst explosions as was suggested by Paczyhski 
(1998), and Woosley (1993). Possible further evidence in support of such a model has 
come from the flattening and reddening of afterglow emission, a few days after the burst, 
in optical wavelength bands (eg. Bloom et al. 1999, Castro- Tirado & Gorosabel 1999, 
Reichart 1999, Galama et al. 1999). 

The goal of this paper is to explore the afterglow flux in different models, uniform ISM 
as well as stratified medium, and compare it with observations in a statistical sense, as 
opposed to comparison with individual GRBs as carried out by Chevalier and Li (1999b). 
We will use this statistical comparison to constrain various physical parameters that 
determine the afterglow luminosity such as the energy, E, the fractional energy in electrons 
(e e ) and magnetic field (e#), the electron energy index (p), and the circumstellar density n. 

In the next section we discuss the afterglow flux and its distribution and comparison 
with observations. 



2. Afterglow flux and its distribution 

Consider an explosion which releases an equivalent of isotropic energy E in a medium 
where the density varies as Ar~ s ; r is the distance from the center of the explosion, and A 
is a constant. The deacceleration radius, r^, where the shell starts to slow down as a result 
of sweeping up the circumstellar material, and the deacceleration time, T^a, in the observer 
frame are given by 
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where T is the initial Lorentz factor of the ejecta, and j3 ~ 1 is a constant. 
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The time dependence of the shell radius and the Lorentz factor can be obtained from 
the self-similar relativistic shock solution given in Blandford & McKee (1976) 
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and t b s is the time in observers' frame at redshift z. 

The magnetic field and the electron thermal Lorentz factor behind the forward shock 
vary as, 
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and 6b and e e are the fractional energies in the magnetic field and the electrons, respectively. 

Using these results we find that the peak of the synchrotron frequency (z/ m ) and the 
cooling frequency (z/ c ), in the observer frame, are 
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The synchrotron self-absorption frequency (in the observer frame) is given by 
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The energy flux at the peak of the spectrum is given by 
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where v v = mm{u m , z/ c } i.e. for u m > u c the peak occurs at v c instead of at v m . 

The equations for s = 2 are as in Chevalier and Li (1999a-b) and are given here for easy 
reference. The flux at an arbitrary observed frequency v can be calculated following Sari, 
Piran and Narayan (1998) in terms of f Up , u m , v c and ua- For the particularly important 
case of v greater than u m and u c the observed flux is: 



f = f u l/2 u (p-l)/2 -p/2 = 
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Note that the flux does not depend on the circumstellar density parameters A and s when 
v > u c , except through an unimportant multiplicative factor (17 — 4s)( p+2 )/ 4 . 

The frequencies v m j '(e^e)/ 2 ), v c t^ and va are shown in figure 1. The afterglow flux for 
the uniform ISM and the wind models only differ when v < v c . Since the cooling frequency 
decreases with time for the uniform ISM model and increases with time for the wind model, 
one of the best ways to distinguish between these models is by observing the behavior of 
the light-curve when v c crosses the observed frequency band as has been pointed out by 
Chevalier & Li (1999). The predictions and comparison of individual GRB lightcurves for 
the two models will be discussed in some detail in a separate paper. Here we turn our 
attention to the statistical property of the afterglow light curve in the two models. 



2.1. Afterglow flux distribution function 

The distribution function for GRB afterglow flux, P t (L u ), at a frequency u, and time 
tobs, g is the probability that the afterglow luminosity (isotropic) is L v at time t b s , g after the 
explosion; z/, t b St9 and L u are measured in the rest frame of the host galaxy. 

The width of Pt(L u ) is a function of the width of the distribution function for E, e e , e#, 
A and p. Assuming that all these variables are independent Gaussian random variables the 
standard deviation (SD) for log(Lj,), ctl v , can be obtained from equation (8), when v > u c 
& u m , and is given by 
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e and E are the mean values of e and E, and oe, &e e , &e B and a p are the standard deviation 
for log-E, loge e , loge# and p respectively; i] for x-ray (10 kev) and optical (2.5 ev) photons 
are shown in figure 2. 

The standard deviation of the flux in the 2-10 kev band at 5 hours after the burst (er^) 
is approximately 0.58 (Kumar & Piran, 1999). Fig 1 shows that this energy band is above 
v c & v m so long as e# > 10 -4 and the density of the surrounding medium is not too small. 
Moreover, rj = 5 for this energy band (see fig. 2), from which we obtain an upper limit 
on a p of 0.26 and the full-with at half-maximum of the distribution for p to be less than 
about 0.6. We note that the electron energy index p lies between 2 and 3 for supernovae 
remanents (cf. Chevalier 1990, Weiler et al. 1986), and Chevalier and Li (1999b) point out 
that the range in p for GRB afterglows is at least ~ 2.1-2.5. 

We can use the variation of 77 with time or v (see figs. 2 & 3) to obtain a p from the 
observed variation to the width of a^. This is of course equivalent to the determination 
of p from the slope of the light curve or the spectral slope. The currently available data, 
consisting of 7 bursts with known redshifts, does not provide a good constraint on the 
temporal variation of a Lv) however, HETE II and Swift would increase the GRB afterglow 
data base by more than an order of magnitude and provide much better constraint on a p . 

Unless €b varies by many orders of magnitude from one burst to another the last term 
in equation (9) is very small and can be neglected. Equating the first two terms individually 
to <j Lv = 0.58 we obtain a E < 0.51, and cr ee < 0.39 (for p = 2.5). For comparison, if the first 
three terms in equation (9) were to contribute equally to a^ v then we obtain Oe = 0.29, 
a €e = 0.22, and o v = 0.15. The mean values for E and e e are not well determined by this 
procedure (but see the discussion below), however E~e e v 1 P+2 can be accurately obtained 
from the observed distribution and is 10 52 erg. 

So far we have discussed a nearly model independent procedure for determining a p and 
a linear combination of a E and of e that relies on making observations in a frequency band 
that lies above u m & u c . In order to determine a E and <r ee separately we need to have some 
knowledge of v c and u m , and therefore the result is model dependent and less certain. For 
instance, if the cooling and the peak synchrotron frequencies are known at some time, even 
if only approximately, then cr^ can be determined from the distribution of the observed 
flux at a frequency v such that v c < v < v m . The SD for fi = L^J^gV 1 / 2 at such an 
intermediate frequency is independent of s and is given by: 



h 16 E 16 eB 16 E 



(11) 
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Once <je is known, equation (9) can be used to determine a ee . Observations at low 
frequencies i.e. v < u m , u c , can be used to constrain a tB and a a which when combined with 
the flux at the peak of the spectrum could be used to determine a tB & a a separately with 
the use of the following equations, 



a l = (4^2^ + I « + O + Q-IYs) °* E > fOT v<Vm< V « ^ 



a l = 9(4 t s) 2 a A + < + ( 12 - 3g ) °b> for v <vc<v m , (13) 



2 1 2 4 2 , /8-3 S \ 2 2 



and 

where /, ^ ^'/'L^ff h - ^V'L,,^- 2 "' 12 - 3 " , f 3 = L„/^'\ and L„, is the 
isotropic luminosity at the peak of the spectrum. 

Another approach to determining the burst and shock parameters is to compare the 
observed flux distributions at several different frequencies and time with the theoretically 
calculated distributions. The latter can be easily calculated by varying E, e e , e#, A and 
p randomly and solving for the flux using the equations given in the last section. Fig. 
3 shows a few cases of flux distribution functions for several different v and t Q b s . The 
advantage of this procedure is that it does not require observational determination of 
various characteristic frequencies i.e. v m , is c , and 1/4, which is difficult to do unless we 
have good spectral and temporal coverage over many orders of magnitude. Moreover, 
this procedure can be used even without the knowledge of burst redshifts; the redshift 
distribution of a representative sub-sample of bursts can be used to calculate the expected 
theoretical distribution of observed afterglow flux which can be directly compared with the 
observed flux distribution to yield the burst and shock parameters. 

It should be noted that if u m , u c , z/4 and L v can be determined accurately then E, e e , 
€b, and A can be obtained for individual bursts as described in eg. Chevalier & Li (1999b), 
and Wijers & Galama (1999), and there is no need to resort to the statistical treatment 
discussed above. 
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3. Conclusion 

We have described how the distribution function for GRB afterglow flux can be used 
to determine the width of the distribution function for the energy in the explosion and the 
shock parameters such as e e , e# (the fractional energy in electrons and the magnetic field), 
p (the power law index for electron energy), and A (the interstellar density parameter). 

The afterglow flux at a frequency above the cooling and the synchrotron peak 
frequencies is independent of interstellar density and scales as E^^^e'f^e^ 2 ^ 4 for 
uniform ISM as well as for energy deposited in a stellar wind with power-law density 
stratification. Using the flux in 2-10 kev band for 7 bursts with known redshifts, 5 hours 
after the burst, which meet the above criteria, we find that the full width at half maximum 
of the distribution for LogE is less than 1.2, for Loge e is less than 0.9 and p is less than 0.6. 
The width for the distribution of p is consistent with the range in p deduced from afterglow 
emissions (cf. Chevalier & Li, 1999b). The width of the distribution for p can be more 
accurately determined from the time variation, or frequency dependence, of the width of 
the afterglow flux distribution. 

For a more accurate determination of the distribution of other parameters we need to 
determine the cooling and the synchrotron peak frequencies (at least approximately), or 
otherwise compare the theoretical and the observed distributions for flux at several different 
frequencies covering the range above and below the cooling and the synchrotron frequencies. 

HETE II and Swift missions are expected to significantly increase the number of 
GRBs with observed afterglow emission and should therefore provide a more accurate 
determination of afterglow luminosity function and the distribution for burst energy and 
shock parameters. 
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Fig. 1. — The figure shows several different frequencies, in electron volts, as a function time 
in observer frame. The solid line is the peak synchrotron frequency for e^ = l&e e = l for 
the wind model i.e. v m / {e l ^ 2 e 2 e ) for s = 0; v m for uniform interstellar medium is larger by 
a constant factor of 1.374. The dotted curve is the cooling frequency, z/ c e^ 2 no, for uniform 
ISM (n is particle number density in the ISM). The dashed curve is u c e % 2 Al for the wind 
model (A* is the baryon density in the wind in the units of 5xlO n g cm -1 ). The energy in 
the explosion E has been taken to be 10 52 erg and p = 2.5. 
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Fig. 2. — r] as defined in equation (10) is shown as a function of time for two different 
frequencies - - 2.5 ev (solid curve) and 10 kev (dashed curve). We took E = 10 52 erg, 
e e = 0.1, and e B = 0.03 for these calculations. 
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Fig. 3. — Distribution function for afterglow flux per unit frequency at photon energy of 1.0 
ev, 0.1 hr after the burst (top left panel), for uniform ISM (solid line) and stratified medium 
with s = 2 (dashed line). The top right panel shows the DFs at photon energy of 100 ev and 
at 0.1 hr after the burst; solid line is for ISM and dashed line is for s = 2. The lower left 
panel shows DFs for photon energy of 1 ev at 1 hr after the burst; solid line is for ISM and 
dashed line is for s = 2. The lower right panel shows the DFs for 10 kev energy photons at 
0.1 hr after the burst (solid line) and at 10 hr (dash-dot curve) for the uniform ISM model; 
the dash-dot curve has been shifted to the right by 2.4. The DFs for the wind model (s = 2) 
for 10 kev energy band are identical to the uniform ISM DFs at times greater than 0.5 hr 
except for an overall shift of the curve to lower energy by a factor of ~ 2. For all of these 
calculations the FWHM of LogE", Loge e , Logee, A and p was taken to be 0.8, 0.7, 3.0, 2.4, 
and 0.25 respectively, and the burst redshift was taken to be 1.5. Moreover, E = 10 52 erg, 
e e = 0.1, Ib = 0.01, and p = 2.3. The diamonds in the lower panel are for bursts with known 
redshifts; the error bar associated with each point is much larger than the size of the symbol. 



